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ABSTRACT 

In the unification scheme of active galactic nuclei (AGN), Seyfert Is and Seyfert 2s are in¬ 
trinsically same, but they are viewed at different angles. However, the Fe Ko emission line 
luminosity of Seyfert Is was found in average to be about twice of that of Seyfert 2s at given 
X-ray continuum luminosity in the previous work dRicci et al.ll2014l) . We construct an accre¬ 
tion disc-corona model, in which a fraction of energy dissipated in the disc is extracted to 
heat the corona above the disc. The radiation transfer equation containing Compton scatter¬ 
ing processes is an integro-differential equation, which is solved numerically for the corona 
with a parallel plane geometry. We find that the specific intensity of X-ray radiation from the 
corona changes little with the viewing angle 9 when 9 is small (nearly face-on), and it is sen¬ 
sitive to 9 if the viewing angle is large (9 > 40°). The radiation from the cold disc, mostly in 
infrared/optical/UV bands, is almost proportional to cos 0 when 9 < 40°, while it decreases 
more rapidly than cos 9 when 9 > 40° because of strong absorption in the corona in this 
case. For seyfert galaxies, the Fe Ka line may probably be emitted from the disc irradiated 
by the X-ray continuum emission. The observed equivalent width (EW) difference between 
Seyfert Is and Seyfert 2s can be reproduced by our model calculations, provided Seyfert Is 
are observed in nearly face-on direction and the average inclination angle of Seyfert 2s ~ 65°. 

Key words: (galaxies:) quasars: general—accretion, accretion discs—black hole physics 


1 INTRODUCTION 


According to the unification scheme of active galactic nuclei 
(AGN), Seyfert Is (Syls) and 2s (Sy2s) are intrinsically same but 
viewed a t different angles , which leads t o diff erent observational 
features l l Antonuccil 1 19931) . iLiu & Wangl fioid) found that the Fe 
Ka line luminosities of Compton-thin Seyfert 2 galaxies are in av¬ 
erage 2.9 times weaker than their Seyfert 1 counterparts. iRicci et al.l 
120141) found that the Fe Ka line luminosities of Syls are about 
twice of those for Sy2s at a given X-ray continuum luminosity 
(10 — 50 keV). The reason is still unclear. One possibility is that 
such difference is caused by anisotropic X-ray emission from these 
sources. Indeed, the orientation dependence of emis sion from AGN 
has been found and studied for a long time. lNemmen & Brothertorl 
d2010h explored the uncertainty of the bolometric corrections of 
quasars with different_viewing angles based on the accretion disc 
models of lHubenv et alj (12000 ). and found that a value of the bolo¬ 
metric luminosity for a quasar viewing at an angle of ss 30° will 
result in ss 30% systematic error, if the emission from the quasar is 
assumed to be isotropic. iRunnoe. Shang. & Brothertonl 120131) ana¬ 
lyzed a sample of radio-loud quasars and found that the quasar lu¬ 
minosity changes with orientation. The sources viewed in the face- 
on direction are brighter than the edge-on sources by a factor of 
2-3. Zh anal d2005i) explained the observed anticorrelation between 
type II fraction and the X-ray luminosity in 2—10 keV based on 


the AGN unification model with only one intrinsic luminosity func¬ 
tion for these two types of AGNs. Recentlv. lDiPompeo et al.| j20l4l) 
used two observed luminosity functions to investigate the intrinsic 
quasar luminosity function with the correction of a simple projec¬ 
tion effect for the anisotropic emission of accretion disc. They con¬ 
cluded that the orientation dependence is the most important one 
among several potential corrections. They also claimed that more 
complex model of anisotropy may strengthen the orientation effect. 


The black hole accretion disc-corona model has been widely 
used to explain both the thermal optical/UV and the power-law 
hard X-ray emission in the spectral energy distributi ons (SEDs) 
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Dullemond 

2001.2002 


iRozanska & Czemvl 12000k iMerloni & Fabianl 
i. The optically thin and geometrically thick hot coronae 
are vertically connected to the both sides of an optically thick and 
geometrically thin accretion disc. It was also suggested that such 
hot coronae may play an important role in launching r elativistic 
jets o bserv e d in X-ray binaries/AGN ( e .g., [M erloni & Fabian 
'2002t ICaol |2004 Izdziarski et ai] l201ll: IWu etalJ 120131: ICac 


2014|). In the accretion disc-corona system, the observed thermal 


optical/UV emission is believed to originate from the blackbody 
radiation of the thin disc passing through the hot corona. A small 
fraction of soft photons are inverse Compton scattered by the 
hot electrons in the corona, which contribute to the observed 
power-law hard X-ray emission of the system. Moreover, the 
temperature of the electron in the transition layer between the 
corona and the disc experiences a rapid decrease from ~ 10 8-9 K 
(the hot corona) to ~ 10 4 ~ 5 K (the cold disc), the thermal X-ray 
line emission may be produced in such transition zone Ixdl2013l ). 
In most previous works, either the cooling rate in the corona or 
the spectra from the accretion flow are calculated assuming the 
corona to be one parallel-plane. The inverse Compton scattering 
is often computed by the Monte Carlo s imulation based on the 
e scape probability method jPozdniakov. Sobol. & Siuniaevlll977l: 
Kawanaka, Kato, & Mineshige! 120081 ; iLiu. Mineshige, & Ohsugal 


2003k ICathoOSh. 


Assuming that a fraction of viscously dissipated energy in the 
disc is transported into the corona (probably by magnetic fields), 
we can calculate the structure of the disc-corona accretion with a 
set of equations of the accretion flow, such as, energy equation, an¬ 
gular momentum equation, continuity equation, and state equation, 
etc., if the model parameters are given. Then, the emitted spectra 
from the accretion flow can be calculated. In this work, we explore 
the angle-dependent emission of the accretion disc-corona system 
in more detail, by solving a set of equations describing disc-corona 
structure and the equation of radiation transfer in the corona si¬ 
multaneously. The change of the spectra with the viewing angle is 
explored in detail. Our results are compared with the X-ray obser¬ 
vations of Syls and Sy2s. The disc-corona model employed in this 
work is briefly described in §2. The calculation method is intro¬ 
duced in §3. We show the results and discussion in §4 and 5. 


2 THE DISC-CORONA ACCRETION MODEL 


of the disc a r = 0.15 is adopted in the calculations. The right side 
of the equation represents the power radiated from the cold disc via 
blackbody radiation, where Tdi sc is the effective temperature in the 
mid-plane of the disc, and r is the optical depth in vertical direction 
of the disc. The continuity equation of the disc is 

- 4tt RH d (R)p(R)vii(R) = M, (3) 


where Hd(R) is the half thickness of the disc, p(R) is the mean 
density of the disc, and dr (R) is the radial velocity of the accretion 
flow at radius R. The state equation of the gas in the disc is 

pkT di 


Ptot — Pgas T Prad — 


/itrtp 


isc -L 4 

' disc} 


(4) 


where p = (1/pi + l/p e ) pi = 1.23 and p e = 1.14 
adopted. The half thickness of the disc Hd is 

y/ptot/p 


Hd — Cs / Hk — 


ftk 


The energy equation of the corona is 

Qcor Q cor 4“ 8Q C or 2 *cor 5 


(5) 


( 6 ) 


where Qcot * s the energy transfer rate from the i ons to the elec¬ 
trons via Coulomb collisions (see equation 11 in ICaol 12009 ) , S 
is the fraction of the energy directly heats the electron, = 
F~ n + F t ~ ern + T’comp is the cooling rate in unit surface area 
of the corona via synchrotron, bremsstrahlung, and Compton 
emissions. The value of 8 can be as high as ~ 0.5 by mag¬ 
netic reconnection, if the mag netic field i n the plasma is strong 
oBisnovatvi-Kogan & Lovelacell 19971 , 120001) . Almost all the power 
dissipated in the hot corona is radiated away locally, which means 
that the radiated power in the corona is independent of the value 5, 
and the temperature and density of the electrons in the corona are 
almost insensitive with this parameter (see lCadl2009l . for the dis¬ 
cussion). In the accretion disc-corona model, the detailed physical 
mechanism for generating the energy source and heating the corona 
is still unclear, though there are some assumptions of corona heat¬ 
ing processes, such as, magneti c fields reco nnection assumed in the 
previous works (e.g. Di Matteolll9 98: Di Ma tteo, Celotti. & Fabianl 
ll999klMerloni & Fabianl2(XriIl2002l:lCad2009!:lXul2()ll) . To avoid 
the complexity, we introduce a parameter, / cor , the ratio of the 
power dissipated in the corona, Q+, r , to the gravitational power 
dissipated in the disc, Qj issi , 


The accretion disc-corona model used in this work is described 
in the previous works dCadl2009l : IXull2013h . The de tailed mode l 
description and calculating approach can be found in lCad 1200 91). 
Here we only briefly summarize the main features of the model. 
The energy equation of the cold thin disc is 

Qdissi - Qtor + \ (1 - ar)Qtor = , (D 

where. 


Q 


+ 

dissi 


^-mqUr) 

07T 



( 2 ) 


is the gravitational power dissipated in unit surface area of the ac¬ 
cretion disc at radius R [where, Qk{R) is the Keplerian velocity, 
Rin = 3f?s, Rs = 2GMdh/c 2 is the Schwarzschild radius for a 
black hole of mass Mbh, and M is the mass accretion rate of the 
black hole]. The third term in the left side of equation (1) represents 
that about half of the power dissipated in the corona is radiated 
back into the disc by Compton scattering, and the reflection albedo 


Qlor = /corQJ issi , (7) 

in our model calculations. 


3 RADIATIVE TRANSFER IN THE CORONA 

We consider a parallel plane geometry of the corona above/below 
the disc. The cooling processes in the corona include Compton, 
bremsstrahlung, and synchrotron emissions. The incident photons 
from the disc at z = Hd are assumed to be blackbody radiation. 
For simplicity, the electron density and temperature of the corona 
are assumed to be constant in the z direction. 


3.1 Radiative transfer equation 


The radiative transfer equation of the corona is, 


dl v (z,p) 

ds 


+ « T ) +ju +jv, 


(8) 
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where ds = dz/p, P = cos 9, 9 is the angle of the photon with 
respect to the vertical direction of the disc, and J„(z, p) is the spe¬ 
cific intensity. For the different absorption and emission processes 
included in the corona, we have the combined absorption and emis¬ 
sion coefficients, = k® + kt, is the total absorption coeffi¬ 
cient including free-free(bremsstrahlung+synchrotron) absorption 
and Thomson scattering, j v = j„ + jf is the emission coeffi¬ 
cient (emissivity) including free-free(bremsstrahlung+synchrotron) 
emission and Comptonization emission. We can calculate Iv{z, p) 
when the structure of the corona, such as electron temperature 
Te,cor(r) and electron density n e , C or(r), are given. 


3.2 Absorption and emission coefficients 

The absorption coefficient of Thomson scattering is kt = 
rie,cor 0 T, where or is Thomson cross section. The absorption co¬ 
efficient including bremsstrahlung and synchrotron processes can 
be described as 

Ku=jv/Bv, (9) 


where B v is the blackbody emissivity, 


B.(T) 


2 hv z /c 2 

exp ( hv/kT ) — 1 


( 10 ) 


and jv = (x® rem + X^ yn )/ 47 r is the emission coefficient of 
bremsstrahlung and synchrotron processes. The bremsstrahlung 
emissivity y^ r em and sync hrotron emissiv ity x^ yn are taken from 
iNaravan & Yil J 1 9951) and lManmotol d2000h . 

The bremsstrahlung emissivity is given by 


brem - ?, ( tlV \ . . 

Xu = %remGexp^ — j, (11) 

where G is the Gaunt factor as in[Rybicki & Lieh tmanl d 19861) . 


„ h (3 kT e \ 1/2 , kT e 

G = W e \nT l -) for ^7 <1 ’ 

h V3 


G = -2- — i n 

kT e 7T 


4 kT e 
( hv 


r kT e - 
for -— > 1 . 
hid 


( 12 ) 

(13) 


The bremsstrahlung cooling rate per unit volume q brem consists of 
electron-ion and electron-electron rates, 


9brem = 9ei + 9ee- (14) 

The electron-ion cooling rate is 

q~i = 1.25nla T ca { m e c 2 F e i(9e), (15) 

where af is the fine-structure constant, 9 e = kT e /m e c 2 is the di¬ 
mensionless electron temperature, and the function F e i has the form 

/ 20 \ 

F ei (6 B ) = 4 (-J-) (1 + 1.T81^e 34 ) for 9 e < 1 , (16) 

Q Q 

Fei(0e) = —^[ln(1.123 6e + 0.48) + 1.5] for 0 e > 1. (17) 

Z7T 

The electron-electron cooling rate is 

q~ e = nlcrlaic (44 - 37 t 2 ) 6> 3/2 (1 + t,19 e + 6> 2 - 1.25#e /2 ) 

for 9 e < 1 , (18) 

q~ = n 2 e cr 2 ai c 2 249 e (ln 1.1232 9 e + 1.28) for 8 e > 1, (19) 

where r e = e 2 /m e c 2 is the classical electron radius. 

The synchrotron emissivity is given by 


= 4.43 x 10 ~ 30 


4Trn e is ,. 47 r m e cv. 

A' 2 (l/ 6 > e ) ( 3 eBdl 


( 20 ) 


and the function I'(x) is given by 


I\x) 


4.0505 0.4 

2-1/6 ( + a-1/4 


— 775 -) exp(—1.8899x 1/3 ), (21) 


where K 2 is the second order modified Bessel function, B is the 
magnetic field strength. For the equipartition magnetic field case, 
we have B = ^/87rp gas ,cor (Pgas,cor is the gas pressure in the 
corona). 

The emissivity of the in verse Compton scattering is calculated 
with the method proposed in lCoppi & Bland ford 119901) . 


■ c 


Ju 


fu(T e , n e ,I v f 
4n 


( 22 ) 


where 


I v i = J I„'(z, f)KT2Trdfi, 


(23) 


is the seed photon intensity of the unit volume including incident 
photons from all the directions. 


3.3 Numerical solution to the radiative transfer equation 

The radiative transfer equation l[ 8 } is an integro-differential equa¬ 
tion. For a given disc-corona structure (including temperatures and 
densities of the electrons and ions, half thicknesses of the disc Ft a 
and thicknesses of the corona H cor , etc.), one can solve this equa¬ 
tion by iterations (e.g.. lCao et al.lfl998h . An initial solution can be 
obtained by solving the equation neglecting the Compton scattering 
term, i.e., 

dFi/(z,/r) _ , w ff , . .fr 

-j- — —Iv{z, p.)(K„ + kt) + ■ (24) 

ds 

We solve the above equation numerically for the corona with the 
boundary condition at the disc surface, 

I u (z = i? d , p) = Bv(Tf), (25) 

in the range 0 < p < 1 , where TJ is the temperature of the disc sur¬ 
face, and the boundary condition at the upper surface of the corona, 

J„(z = Hd + Ft cor , ft) = 0, (26) 


in the range — 1 < p < 0. With derived initial solution, we cal¬ 
culate the emissivity of the inverse Compton scattering from equa¬ 
tions < 1221 and 031 . and solve equation { 8 } numerically. With de¬ 
rived I v (z, p), the emissivity of the Compton emission can be re¬ 
calculated with equations 1221 and ( 123b . We find that the final so¬ 
lution can be achieved after several iterations when the solutions 
converge. 

Integrating the intensities over the different directions (p = 
cos 9) and frequencies u, the energy loss via radiation from the 
upper and lower surfaces of the corona can be calculated by 


Qcor = du 


Iv(z = Hd + H COI , n)2npdp 


+ J dv J I v (z = Hd, p)2-k pdp. (27) 


Subtracting the incident blackbody radiation from the thin disc, 


<2disc= / KBvdv, 


(28) 
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we obtain the cooling rate in unit surface area of the corona, 

F~ t = QSr - Qa isc ■ (29) 


Given the black hole mass Mbh, the dimensionless mass ac¬ 
cretion rate rh (rh = M / M^dd, where Mvdd = TEdd/O.lc 2 ), 
and the fraction of the energy directly heats the electron 8, we can 
derive the disc structure (such as, the effective temperature in the 
mid-plane of the disc, the density in the disc, and the half thick¬ 
ness of the disc, etc.) as a function of radius R from equations 
0 -©. and equation 0 . when the ratio of the power dissipated 
in the corona / cor is specified. The structure of the corona (such as, 
temperatures and densities of the electrons and ions, scaleheight 
of the corona) can be derived with equations © and 1291 under 
the assumption of equipartition of the magnetic pressure and the 
gas pressure in the corona. We assume the temperature of the ions 
in the corona T i cor — 0.9T v i r = 0.9GMm p /3kR in this work 
as that in lCaol 120091 ). Thus, the accretion disc-corona structure is 
available by solving the radiation transfer equation together with 
disc-corona equations described in § 2 . 

The specific intensity from the corona at a certain radius R is 
obtained as a function of direction and frequency of the photons. 
Integrating the intensity over the whole surface of the corona, the 
specific luminosity emitted from the corona per steradian is 


dLii 
2n d/j, 


= / iih(z = H d + .Hcor, /r) 27 rrdr. 


(30) 


Thus, we can obtain the direction-dependent spectrum of an accre¬ 
tion disc-corona system for a set of given disc parameters. 


4 RESULTS 

We adopt the model parameters as follows, the black hole mass 
Mbh = 10 s M 0 , the maximum radius of the disc-corona flow 
f?cormax = 100f?s, and the fraction of energy directly heated the 
electron 5 = 0.5, are fixed in all cases. The mass accretion rate rh 
and the ratio of the power dissipated in the corona / cor are taken as 
two free parameters in our model calculations. 

In Figure 1, we show the emergent spectra for four disc-corona 
accretion models with different values of parameters rh and / cor . 
The values of the two parameters in the models are, model (a): 
rh = 0.1 and / cor = 0.1, model (b): rh = 0.1 and / cor = 0.3, 
model (c): rh = 0.5 and / cor = 0.06, and model (d): rh = 0.5 
and / CO r = 0.1. For each model, we plot the emergent spectra ob¬ 
served in five different directions with the lines in different types 
and colors. The spectra from 6 = 87°, 84°, 76°, 50°, and 2° are 
shown with black dashed, cyan dotted, magenta dash-dotted, green 
thin solid, and blue thin dash-dotted lines, respectively. The corre¬ 
sponding value of 8 is denoted near each spectrum. The red thick 
solid line is the spectra integrated over the directions from /j, = 0 
to /r = 1. The four vertical dotted lines represent the four typical 
frequency points corresponding to 2500A, 0.1 keV, 2 keV, and 10 
keV, respectively. We can see from these figures that the spectra 
from different emitting directions are different both in the luminos¬ 
ity and the spectral shape. The luminosity decreases with the in¬ 
creasing viewing angle between the emitting photons and the axis 
of the accretion disc, 8. The spectral shapes are similar in almost 
all bands but very different between the optical/UV (~ 2500A) and 
soft X-ray band (0.1 ~ 1 keV). 

The different values of mass accretion rate rh = 0.1 and rh = 
0.5 are adopted in Figure 1(a) and Figure 1(d), respectively, while 
the values of all other parameters are the same in these two figures. 


The larger the mass accretion rate adopted, the stronger and harder 
the spectra are. The different values of / cor = 0.1 and / cor = 0.3 
are adopted in Figure 1(a) and Figure 1(b), respectively, while the 
values of all other parameters are the same in these two figures. We 
find that the change of the spectra with / cor is also evident. The 
larger the value of / cor adopted, the harder the spectra are. 

The spectral shapes plotted in the four panels of FigureQ]show 
a certain degree of degeneracy of the two parameters rh and / cor . 
Different combinations of these two parameters may give very sim¬ 
ilar spectral shapes, but these can be distinguished by the different 
luminosity, which is directly controlled by rh. 

In order to study the change of spectra with the viewing direc¬ 
tion, we calculate some typical quantities of the observed spectra, 
including the observed bolometric luminosity Lboi, the spectral lu¬ 
minosity at optical/UV band (A=2500 A) L 0 , and the spectral lu¬ 
minosity at X-ray band (E — 2 keV) L x . The optical/UV to X-ray 
power index a ox is defined as 

1OgL ,(2 5 00 

log ^(2500 A) /v(2 keV) ’ 

and the X-ray spectral index between 2 keV and 10 keV, a x (de¬ 
fined as iJ^x 10 keV ' ) oc t/ _ “ x ). The changes of these quantities 
with the viewing angles are plotted in Figures 2-5. 

In Figure 2, we plot the change of observed bolometric lumi¬ 
nosity with viewing angle. The red solid, green dashed, cyan dotted, 
and black dash-dotted lines correspond to the four models (a)-(d), 
respectively. The blue thin dashed line indicates a simple relation 
as Lbol oc cos 8 representing the area-projection effect. 

In Figures 3 and 4, we plot the changes of observed spectral 
luminosities at optical/UV band (A = 2500A) and X-ray (2 KeV) 
band with the viewing angle. The red solid, green dashed, cyan dot¬ 
ted, and black dash-dotted lines correspond to the four models (a)- 
(d), respectively. The blue thin dashed lines also represent simple 
relations as L a oc cos 8 and L x oc cos 8. We find that the two spec¬ 
tral luminosities (optical and X-ray bands) show different changes 
with the viewing angle 8. The X-ray luminosity L x decreases more 
slowly than L x oc cos#, while L a decreases more rapidly than 
L 0 oc cos 8. The X-ray luminosity L x is almost isotropic when the 
viewing angle is small (nearly face-on), and it becomes strongly 
anisotropic if the viewing angle is large (> 30° — 40°). The optical 
luminosity L a is almost proportional to cos 8 when 8 < 30° — 40°, 
but decreases more rapidly than cos 8 when it is viewed at large 
angles (> 30° - 40°). 

The changes of the optical/UV to X-ray spectral index, a ox , 
and the X-ray spectral index between 2 keV and 10 keV, a x , with 
the viewing angle are shown in Figure 5. The red solid, green 
dashed, cyan dotted, and black dash-dotted lines correspond to the 
four models (a)-(d), respectively. In each model, a ox decreases very 
slowly with 8 at small viewing angles and quickly at large angles 
nearly edge-on. The situation is different for a x , which is almost 
unchanged with 8 in each model. 

The hard X-ray continuum emission is from the coronae above 
the discs, and the discs are irradiated by the X-ray photons from 
the corona. Our model calculations show that the angle-dependent 
X-ray continuum spectra are anisotropic. They deviate from cos in¬ 
dependence. The Fe Ka lines are probably emitted from the irra¬ 
diated discs. In this case, the Fe line emission is anisotropic, and 
its angle-dependence follows ~ cos 8. Thus, we can calculate the 
observed equivalent widths as functions of viewing angle with dif¬ 
ferent values of model parameters. The result of model (a) is shown 
with red solid line in Figure 6 . 
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It is still unclear whether part of observed Fe line emission is 
from the torus, which is emitted nearly isotropically. We estimate 
how the results would be affected by the toms contribution to Fe 
line emission by assuming that x per cent of the total Fe line lumi¬ 
nosity is from the toms. The ratio of the equivalent widths of the 
Fe line emission viewing at an angle 6 is 

EW{9) _ (100 — X) COS0 + X ^ L v (2 keV)(0°) 

EW(0°) ~ 100 X L„(2k„v)(0) ' ( } 

We re-calculate the relative equivalent width of the Fe line for dif¬ 
ferent values of a;=50, 20, and 10. The results are shown in Figure 
6 with black dashed, green dotted, and blue dash-dotted lines, re¬ 
spectively. 


5 DISCUSSION 

Solving the radiative transfer equation of the corona numerically, 
we obtain the cooling rate in the hot corona and the emergent spec¬ 
trum of the accretion disc-corona system viewed at any specified 
angle. The calculations are carried out for four sets of model pa¬ 
rameters, in which different accretion rate rh and different ratio of 
the power dissipated in the corona / cor are adopted, but all the other 
model parameters are fixed. Overall, the calculated spectra are sen¬ 
sitive to the values of the model parameters. For the models with 
same value of / cor [compare Fig. 1(a) with 1(d)], the larger the ac¬ 
cretion rate, the harder the spectra from the disc-corona accretion 
flow are. The X-ray emission predominantly originates from the in¬ 
verse Compton scattering of the soft photons from the thin disc by 
the hot electrons in the corona. The larger the accretion rate, the 
more the soft photons are radiated from the disc, thus the emission 
in the lower energy band (i.e. optical/UV band) becomes relatively 
weaker while the high energy band (i.e. X-ray band) becomes rela¬ 
tively stronger, and the spectra are harder (see a ox in the left panel 
of Fig. 5). On the other hand, for the models with same value of rh 
[compare Fig. 1(a) with 1(b), or Fig. 1(c) with 1(d)], the higher 
the ratio of the power dissipated in the corona, / cor , the harder 
the spectra from the disc-corona accretion flow are. We find that 
the electron temperature of the corona increases with / cor , which 
causes the Compton emission to be hard for the high / cor cases. 

It is found that the spectra of disc-corona system observed at 
different directions have the similar shapes in 2-10 keV X-ray band. 
The X-ray spectral index, a x , remains almost unchanged for differ¬ 
ent viewing angle 9 (see the lower panel of Figure 5). The X-ray 
emission mostly originates from the inverse Compton scattering of 
the soft photons radiated from the thin disc by the hot electrons in 
the corona. The inverse Compton scattered X-ray spectrum mainly 
depends on spectrum of the seed photons and the temperature of the 
hot electrons. Thus, the X-ray spectral index is not sensitive to the 
viewing direction. Comparing the values of a x for different models, 
we find that a x « 1.64 for the case with rh = 0.1 and / cor = 0.1, 
while a x ~ 1.40 for the case with rh = 0.5 and / cor = 0.1. This 
result seems inconsistent with the observed results and theoretical 
predictions in the previous works which claim the increasing a x 
with the Eddington ratio, Lboi/TEdd or accretion rate, rh (e.g., see 
Fig. 4 in iCacll2009l ). The reason is that we employ the same value 
of parameter, / cor = 0 . 1 , for these two models, which is inconsis¬ 
tent with the fact that the value of / cor decreases with the accretion 
rate (e.g., see Fig. 1 in ICaol 12009 ). If we compare the results of 
the models (b) with (c), which adopt / cor = 0.3 for the case with 
m = 0.1 and a smaller f COI = 0.06 for the case with rh = 0.5, we 
have q x ~ 1.43 for the case with rh = 0.1 and / cor = 0.3 (see 
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Figure 1. Spectra of the four disc-corona accretion models with different 
values of parameters rh and / CO r in four plots, model (a): rh = 0.1 and 
/cor = 0.1, model (b): rh = 0.1 and / cor = 0.3, model (c): rh = 0.5 and 
f cor = 0.06, and model (d): rh = 0.5 and f co r = 0.1. For each model, 
we plot the spectra observed at different angles (the dashed lines in different 
colors). The emergent spectra viewed at 6 = 87°, 84°, 76°, 50°, and 2°, 
are shown with black dashed, cyan dotted, magenta dash-dotted, green thin 
solid, and blue thin dash-dotted lines, respectively. The corresponding value 
of 6 is denoted near each spectrum. The red thick solid lines are the spectra 
integrated over all directions (p = 0 — 1). The four vertical dotted lines 
represent the four typical frequencies corresponding to 2500A, 0.1 keV, 2 
keV, and 10 keV, respectively. 
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Figure 2. The observed bolometric luminosity varies with the viewing an¬ 
gle. The red solid, green dashed, cyan dotted, and black dash-dotted lines 
correspond to the four models (a)-(d) respectively. The blue thin dashed line 
indicates the relation of Lbol ' x cos 9. 



Figure 3. The spectral luminosity at optical/UV band(A = 2500A) varies 
with the viewing angle. The red solid, green dashed, cyan dotted, and black 
dash-dotted lines correspond to the four models (a)-(d) respectively. The 
blue thin dashed line represents the relation of L a oc cos 9. 


green line Figure 5), while a x ~ 1.62 for the case with rh = 0.5 
and / C or = 0.06 (see cyan line in Figure 5). 

The main focus of this work is to explore the change of the 
emergent spectra from the disc-corona system viewed at different 
angles. In Figure 1, we find that the luminosity decreases with the 
increasing viewing angle 0 with respect to the axis of the accretion 
disc. The detailed results of the observed bolometric luminosity are 
plotted in Figure 2, which shows that the bolometric luminosity is 
nearly proportional to p (p=cos0). This is due to the area-projection 
effect. The change of the spectral shapes for different viewing an¬ 
gles is different between optical/UV and soft X-ray bands (see Fig¬ 
ure 1). The observed spectral luminosity at a typical optical/UV 
band of 2500A decreases with 9 (see Figure 3). It decreases more 
rapidly than cos ^-relation. The blackbody emission from the thin 
disc is partly absorbed in the corona. For the optical/UV emission 
from the disc, the specific intensity of the photons from the upper 
surface of the corona, I„ ~ 1^,0 exp(—ro//x), where o is the 
specific intensity of the photons injected in the lower surface of the 
corona, and to is the vertical optical depth of the corona. Thus, the 
observed spectral luminosity at optical/UV band is proportion to 
^texp(—To/p,), which decreases more quickly with 9 than cos0- 
relation. On the other hand, although the X-ray spectral shape re¬ 
mains unchanged for different 9, the observed spectral luminosity 


Figure 4. The spectral luminosity at X-ray band) E = 2 keV) varies with 
the viewing angle. The red solid, green dashed, cyan dotted, and black dash- 
dotted lines correspond to the four models (a)-(d) respectively. The blue thin 
dashed line represents the relation of L x oc cos 9. 



0 (degree) 

Figure 5. The optical/UV to X-ray spectral index, a ox (upper panel), and 
the X-ray spectral index between 2 keV and 10 keV, « x (lower panel), vary 
with the viewing angle. The red solid, green dashed, cyan dotted, and black 
dash-dotted lines correspond to the four models (a)-(d) respectively. 

at 2 keV decreases with 9 (see Figure 4). It decreases more slowly 
than cos 0-relation. 

The narrow 6.4 KeV Fe Ka lines are ubiquity in AGNs. Al¬ 
though its origin is still uncertain, it was suggested that they may 
probably originate from the distant molecular cloud (torus), the 
outer accretion disc or/and the broad-line region (BLR). The BLR 
origin is ruled out by the fact that no correlation between the Fe Ka 
core width and the BLR line (i.e., H/3) width jNandrall2006t) . The 
emitted narrow 6.4 KeV line are found to show different properties 
for different types of AGN. The line luminosities of the narrow 6.4 
KeV Fe Ka line from type I AGNs are much stronger than those 
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0(degree) 

Figure 6. The equivalent width of the narrow Fe Ka line emitted from 
AGN varies with the viewing angle, which is calculated with model (a). 
The red solid line represents the case that all the line emission is radiated 
from the thin accretion disc. The black dashed, green dotted, and blue dash- 
dotted lines are the results of the cases that 50%, 20%, and 10% of the 
line luminosity contributed from the torus, respectively. The black dotted 
horizontal line represents the relative equivalent width being 0.5. 


from type II AGNs at the same X-ray continuum luminosity. Com¬ 
piling 89 Seyfert galaxies and using [O I V] emission to e stimate 
the intrinsic luminosity of the sources, lliu & Wang ( 2010 ) found 
that the Fe Ka line luminosities of Compton-thin Seyfert 2 galaxies 
are in avera ge 2.9 times weaker than their Seyfert 1 counterparts. 
iRicci et aD ' 2014 b found that the Fe Ka line luminosity is corre¬ 
lated with the 10-50 keV X-ray continuum luminosity either for 
Syls and Sy2s. The slopes of the correlations are almost same for 
these two types of the sources, but the Fe Ka line luminosities of 
Syls are about twice of those for Sy2s at a given X-ray continuum 
luminosity. 

It is believed that type 1 Seyfert galaxies are intrinsically sam e 
as type 2 Seyferts but viewed at different angles j Antonuccl 19931) . 
We find that the observed systematical difference of EW of Fe Ka 
emission lines between Syls and Sy2s can be attributed to the dif¬ 
ference of the view angles(see Figure 6). Such EW difference can 
be reproduced by our model calculations, provided Syls are ob¬ 
served in nearly face-on direction and the average inclination angle 
of Sy2s ~ 65°, which support the unification scheme of AGN. If 
a fraction of Fe line emission is contributed by the torus, a larger 
average inclination angle is required for Sy2s, which implies that 
the contribution from the torus should be much less than that from 
the disc. 

No significant correlation is found between the spec¬ 
tral index a ox and the radio core dominance parameter R, 
which is_ believed to be an indicator of the viewing angle 
feunnoe, Shane, & Brothertonl l2013h . Our results show that the 
spectral index a ox is weakly dependent on the inclination angle 
(see Figure 5). This is not surprising, because there is a strong cor- 
relation between q ox and optical luminosity/Eddington ratio (e.g.. 


Vignali, Bra ndt. & Sch neider 20031: iGrupe et al.ll20ld : lLusso et alJ 


201 A . This correlation may smear out any possible correlation be¬ 

tween a ox and inclination angle for a normal AGN sample. We 
suggest that the investigation on a ox and inclination angle relation 
should be carried out with a sample of AGNs in the narrow range 
of Eddington ratio/luminosity. 

In the standard unification model of AGNs, the different ob¬ 
servation properties of different types of AGN, such as, different 
continuum spectral shape and broad emission lines, can be ex¬ 
plained by their different inclination angles of the black hole ac¬ 


cretion flow and the surrounding torus. We show in this work that 
the emission spectra emitted from the accretion flow of a black hole 
are anisotropic along the different viewing angles, which is consis¬ 
tent with the area-projection effect in the bolometric luminosity, but 
very different in the other characteristics of SED, i.e., the spectral 
luminosities in the optical/UV band and the X-ray band. Our calcu¬ 
lations of accretion disc-corona spectra may provide more precise 
orientation effect correction in predicting the intrinsic luminosity 
functions of AGN sources from the observed luminosities at certain 
band. Recently, the effect of anistropic radiation from the accretion 
discs on the luminosity function derived from an AGN sample was 
evaluated bv lDiPompeo et all d2014l) . and they found that the bright 
end of the luminosity function may be overestimated by a factor 
of ~ 2 without considering this effect. A simple cos ((-dependent 
specific intensity from a bare accretion disc is used in their esti¬ 
mates. Our present detailed calculations of the disc corona spectra 
as functions of the viewing angle can be incorporated in deriving 
the intrinsic AGN luminosity. 

The calculations of the radiation transfer in the corona in this 
work are carried out in Newtonian frame. A general relativistic ac¬ 
cretion corona model is required for direct modeling the observed 
spectra of AGN. The calculations in this work can be easily ex¬ 
panded for the accretion discs surrounding Kerr black holes in gen¬ 
eral relativistic frame. This will be reported in our future work. 
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